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Observation of a narrow baryonic state in DIS at HERA 
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A resonance search has been made in the Kg p, K% p and K'^p invariant-mass spectra measured 
with the ZEUS detector at HERA using an integrated luminosity of 121 pb~^. The search was 
performed in the central rapidity region of inclusive deep inelastic scattering at an ep centre-of-mass 
energy of 300-318 GeV for exchanged photon virtuality, Q^, above 1 GeV^. The results support 
the existence of a narrow state in Kg p and Kg p decay channels, consistent with the pentaquark 
prediction. No signal was found in the K^p decay channel. 
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I. INTRODUCTION 



Recent results from fixed-target experiments triggered 
new interest in baryon spectroscopy after observing a 
narrow baryonic resonance in the K^n decay channel 
with a mass of approximately 1530 MeV and positive 
strangeness 0, S S 13 ■ Such a state can be interpreted 
as a bound state of five quarks, i.e. as a pentaquark, 
0+ ~ uudds. According to the predictions of the chiral 
soliton model ^,K%p and K% p (denoted as K° p (p)) 
decays are also possible. For the Kg p channel, evidence 
for a corresponding signal has been seen by other low- 
energy fixed-target experiments 0, H S IS ^ • Recently, 
two other pentaquark candidates have also been reported 

urn 

The Kg p [p] decay channel is expected for the PDG 
E bumps 13], which have been observed by early fixed- 
target experiments. Such unestablished resonances are 
too low in mass to be accommodated in most quark mod- 
els. The presence of such states in the mass regions close 
to the production threshold of Kg p (p) final state com- 
plicates the search for possible pentaquark. 

The 6+ state and the S bumps discussed above have 
never been observed in high-energy experiments, where 
hadron production is dominated by fragmentation. 

This paper discusses recent ZEUS results 14] on the 
search for exotic baryons in the Kg p (p), K^p, K^p 
invariant-mass spectra measured with the ZEUS detector 
at HERA. The measurements were based on the central 
pseudorapidity region {\ri\ < 1.5) where the contribution 
from the fragmentation of the proton remnant is negligi- 
ble. If the pentaquark state is produced without the net 
baryon number carried by the proton remnant, emerging 
from the emissions of giuons and quarks in the hadro- 
nisation process, this would open a new chapter in our 
understanding of non-perturbative QCD. 
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FIG. 1: The tt^tt" invariant-mass distribution for > 
1 GeV^. The solid line shows the fit result using two Gaus- 
sians plus a linear background, while the dashed line shows 
the linear background. 



II. EVIDENCE FOR A BARYONIC STATE 
DECAYING TO A'g p (p) 

The present analysis was performed using deep inelas- 
tic scattering events measured with exchanged-photon 
virtuality > 1 GeV^. The data sample corresponds 
to an integrated luminosity of 121 pb~^, taken between 
1996 and 2000. The analysis is based on charged tracks 
measured in the central tracking (CTD). The tracks were 
selected with pt > 0.15 GeV and \ri\ < 1.75, restricting 
this study to a region where the CTD track acceptance 
and resolution are high. 

X^-mesons were identified by their charged-decay 
mode. Kg tt+tf" (see Fig. The candidates with 
PT(ifs) > 0-3 GeV and l?7(ifg)l <_1.5 were retained. To 
eliminate contamination from A(A) decays, candidates 
with a proton mass hypothesis AI{ptt) < 1121 MeV were 
rejected. 

The (anti)proton-candidate selection used the energy- 
loss measurement in the CTD, dE/dx. As example. 
Fig. 121 shows the dE/dx distribution as a function of the 
track momentum for negative tracks. Clear antiproton 
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FIG. 2: The dE/dx distribution as a function of the track 
momentum, p, for negative tracks. The sohd hues indicate 
the antiproton band used in this analysis, while the dotted 
lines indicate the selection cuts dE/dx > 1.15 mips and p < 
1.5 GeV used in the Kg p (p) reconstruction. 



and iC" -meson bands were observed (dE/dx for Monte 
Carlo simulations and for positive tracks have similar 
qualities, not shown). 

Kg p (p) invariant masses were obtained by combin- 
ing Kg candidates in the mass region 480 — 510 MeV 
with (anti)proton candidates in the (anti)proton dE/dx 
band with the additional requirements p < 1.5 GeV and 
dE/dx > 1.15 mips in order to reduce the pion back- 
ground. The CTD resolution for the Kg p (p) invariant- 
mass near 1530 McV, estimated using Monte Carlo sim- 
ulations, was 2.0 ± 0.5 MeV for both the Kg p and the 
Kg p channels. 

Figure |21 shows the Kg p (p) invariant mass for > 
20 GeV^ , as well as for the Kg p and Kg p samples sepa- 
rately (shown as inset). The distribution of the Ariadne 
Monte Carlo model was normalised to the data in the 
mass region above 1650 MeV. The data are above the 
Monte Carlo prediction near 1470 MeV and 1522 MeV, 
with a clear peak at 1522 MeV. 

To extract the signal seen at 1522 MeV, the fit was 
performed using a background function plus two Gaus- 
sians. The background has the threshold form F{M) = 
Pi{M — rriK — rup)^^ (1 -t- P^{M — rriK — i^p)): where rriK 
and nip are the masses of the kaon and the proton, respec- 
tively, and /i=i,2,3 are free parameters. The first Gaus- 
sian, which significantly improves the fit at low masses, 
may correspond to the unestablished PDG S(1480). 
The second Gaussian describes a pronounced peak at 
1522 MeV. The peak position determined from the fit 
was 1521.5 ± 1.5(stat.)t^;?(syst.) MeV. It agrees well 
with the measurements by HERMES, SVD and COSY- 
TOF for the same decay channel 0, IS ■ The measured 
Gaussian width of 6.1 ± 1.6(stat.)t2^4(syst.) MeV was 
above, but rather close to the experimental resolution. 

In order to determine the natural width of the pos- 
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FIG. 3: Invariant-mass spectrum for the K% p (p) channel 
for > 20 GeV^. The solid line is the result of a fit to 
the data using the threshold background plus two Gaussians. 
The dashed lines show the Gaussian components, while the 
dotted line indicates background. The prediction of the Ari- 
adne simulation is normalised to the data in the mass region 
above 1650 MeV. The inset shows the Kgp (open circles) and 
the Kg p (black dots) candidates separately, compared to the 
result of the fit to the combined sample scaled by a factor of 
0.5. 
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Gaussian -l-Bkg. 


1522.2 ± 1.5 


4.9 ± 1.3 


155 ± 40 


51/47 


2 Gaussians-|-Bkg. 


1521.5 ± 1.5 


6.1 ± 1.6 


221 ±48 


35/44 



TABLE L Fit resuhs (in MeV) for > 20 GeV^ All 
parameters are free in the fitting procedure. 



sible state, the Breit-Wigner function convoluted with 
the Gaussian was used in the fitting procedure. If the 
width of the Gaussian is fixed to the experimental res- 
olution, the estimated intrinsic width of the signal was 
r = 8 ± 4(stat.) MeV. The systematical error on this 
width is expected to be smaller than the statistical error, 
but due to low statistics, complicated background and 
the narrowness of the peak leading to unstable fits, full 
systematical uncertainty was difficult to estimate. 

The number of events ascribed to the signal by this 
fit was 221 ± 48. The statistical significance, estimated 
from the number of events assigned to the signal by the 
fit, was 4.6cr. Other fit parameters are shown in Table 1. 
The number of events in the Kgp channel was 96 ± 34. If 
the signal corresponds to the pentaquark, this provides 
the first evidence for its antiparticle. 

Under a conservative assumption, one could assume 
that the enhancement from the left side of the 1522 MeV 
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FIG. 4: Same data as in Fig. |3 but fit is performed using a 
single Gaussian plus the background. The x^/ndf is accept- 
able since it is dominated by the region > 1550 MeV. For 
M < 1550 MeV, the xVpoint is 32/22 for the same fit. 



peak has a statistical nature, or is due to an instrumental 
effect which is not taken into account in the Monte Carlo 
simulation. As seen in Fig.^ the single-gaussian fit plus 
background is still acceptable according to the x^/ndf 
test. However, in the low-mass region, M < 1550 MeV, 
this fit had a poor quality. For the background fit only, 
the x^/point = 2.2, while the overall x^/ndf, which is 
dominated by the high-mass region, is reasonable (see 
Table 1). 

It is interesting that the K'^ p data has a peak near 
1470 MeV, while the Kg p channel does not have a clear 
peak at this mass. Generally, one could conclude that the 
data may contain a contribution from the unestablished 
I](1480) state. This conclusion is mainly supported by 
the comparison of the data with the Monte Carlo simu- 
lation. The parameters and the significance of this state 
are difficult to estimate due to the steeply falling back- 
ground. 

An ensemble of Monte Carlo experiments was gener- 
ated in which the background shape was parameterised 
by the threshold function, and the background rate was 
simulated as the mean of a Poisson distribution in each 
mass bin. The probability of a similar signal anywhere in 
the range 1500-1560 MeV arising from statistical fiuctu- 
ation of the background was below 6 x 10~^. For a more 
realistic case, when the background was parameterised 
by the threshold function plus the 1465 MeV Gaussian 
as the starting distribution, the probability was found to 
be about a factor of ten lower. 

The mass spectrum was studied over a large range. 
The signal was well seen for > 60 GeV^, but, in 
this case, the statistics is smaller. No significant signal 
was found at very low {Q^ = 1 — 5 GeV^) due to a 
complicated background and possible acceptance effects. 
However, the same signal was found at > 1 GeV^ at 
low photon-proton centre-of-mass energy, W < 120 GeV. 
For this W range, the charged multiplicity (and, hence. 



FIG. 5: Invariant-mass spectra for the K'^p and K p chan- 
nels (plus charge conjug ates) for > 1 GeV^ 



the combinatorial background) is lower. 

In conclusion, the results constitute evidence for the 
existence a narrow baryonic resonance decaying to the 
Kg p (p) channel with the mass near 1522 MeV. This is 
the first observation of such resonance in high-energy col- 
liding experiments in the phase-space regions dominated 
by fragmentation. 

The measured mass and the width of the observed 
state are close to those observed in the K~^n channel Tij, 
and agree well with the theoretical expectations 5] . The 
PDG reports no any E state in the invariant-mass region 
1480-1560 MeV, and no peak with a similar mass was 
observed in Att decays. This favours the conclusion that 
the observed peak can be interpreted as a pentaquark 
candidate. In this interpretation, the signal seen in the 
Kg p channel corresponds to an antipentaquark with a 
quark content of uudds. 



III. e++ STATE? 

If 8+ is an isotensor state, a 0++ signal can be ex- 
pected in the K'^p spectrum The K^p{K^p) in- 
variant mass spectra were investigated in a wide range of 
minimum values, identifying proton and charged kaon 
candidates using the dE/dx information. The proton 
candidates inside the dE/dx proton band were required 
to have dE/dx > 1.8 mips, while the kaon candidates 
were reconstructed in the kaon band after the restriction 
dE/dx > 1.2 mips. For > 1 GeV^, no peak was 
observed near 1522 MeV in the K~^p and K~p spectra, 
while a clean signal was seen in the K~p{K^p) channel 
at 1518.5 ± 0.6(stat.) MeV, corresponding to the PDG 
A(1520)Do3 state, Fig.gl 

The numbers of the reconstructed A(1520) and 
A(1520) candidates agree very well with each other. The 
same good agreement was found for K^ and K~ , as well 
as for p and p candidates. This indicates negligible contri- 
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FIG. 6: Invariant-mass spectra for the K^p and K~p chan- 
nels for Q2 > 20 GeV^ 



butions from the proton remnant, as well as from possible 
secondary-scattering events. 



The K'^p and K~p mass spectra were investigated at 
> 20 GeV^, i.e. in the kinematic region where the 
1522 MeV peak is clearly pronounced for the Kg p (p) 
channel. Fig. El shows the mass spectra for > 
20 GeV^. Again, no sign of a peak was found in the 
K^p decays. 



The failure to observe 0+^ 
is not isovector or isotensor. 



indicates that the 9+ state 



In contrast to the Kg p (p) state, the detector ac- 
ceptance for A(1520) baryons decreases faster with in- 
crease of , since the momenta of ii'^-mesons were re- 
stricted by the requirement dE/dx > 1.2 mips (i.e. high- 
momentum mesons were not used in the reconstruction), 
while K'g used in the K'g p (p) reconstruction were not 
restricted in the total momenta. Since the detector ac- 
ceptances are significantly different for the K^ p (p) and 
K~p measurements, direct comparisons of the produc- 
tion rates of possible pentaquarks and A(1520) state are 
not possible before taking into account the detector ac- 
ceptance effects. 
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